Strong Coupling Theory of the Spinless Charges on the Triangular Lattices: 

Possibility of a New Quantum Liquid 
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We propose a new type of charge liquid state in the spinless fermion system on a triangular lattice 
under strong inter-site Coulomb interactions, V. In the strong coupling limit (t — 0), the ground 
state is classical and disordered due to geometrical frustration. The introduction of small t will 
drive the system to a partially ordered phase which we call a "pinball liquid". A possibly long range 
ordered Wigner crystal solid coexist with a liquid component which are moving around them like a 
pinball. This liquid is dominant over wide range of filling, even away from the regular triangle, and 
is also realized in the hard core boson systems. Relevance to the organic 0-ET2X is discsussed. 

PACS numbers: 75.40, 75.50, 72.80.Sg, 72.80.Le 
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Geometrical frustration induces exotic states in the 
strongly correlated systems due to the severe competi- 
tion among ordered states. They are explored exper- 
imentally in many materials from transition metals like 
NiGa 2 S 4 [3, Na^CoOaH and YFe 2 4 |3 to organic solids, 
k-ET 2 Cu 2 (CN) 3 Q. Theoretical studies on frustrated 
spin systems arc conventionally intented on a search for 
" spin liquids" [f| which generally lacks any long range or- 
der. However, although the magnetic order is destroyed 
by frustration another non-magnetic type of order often 
appears instead. Typical example is a dimer formation^ 
and recently, a new gapless state called nematic order is 
found in the vicinity of ferromagnetic phase |7],|8j. 

Charge degrees of freedom under frustration has also 
been an important topic from the past0, and one of the 
latest topics is a melting of charge order(CO) into the 
metallic state on the triangular lattice^3| m analogy with 
the "spin liquid" . Although this seems to support an "or- 
derless charge liquid" picture, the possibility of another 
type of ordering still remains unclear. In this Letter, 
we propose a new type of partially ordered charge liquid 
in such triangular lattice system. This liquid is nobel 
in that a charge degrees of freedom shows a coexistence 
of "statically ordered" and "liquid" components, where 
a gapless spatial long range order is possibly realized. 
The results are to be compared with an inhomogeneous 
metallic state found in 6»-ET 2 CsZn(SCN) 4 [ll|. 

One of the typical cases that CO resides is the 
quarter-filled electronic systems under strong electronic 
interactions 0| . It is conventionally described by 
the extended Hubbard model(EHM) including the on- 
site and inter-site Coulomb interactions, U and V, 
respectively^ besides the transfer integrals, t. Since 
the spin degrees of freedom plays only a secondary role 
when t < V < U, i.e. J ~ t 2 /U <C t, treating the 
spinless fermions at half-filling, p=0.5, is enough to un- 
derstand the CO phenomena realized in the temperature 
range of J < T < Too, where Tco denotes the onset 
temperature of CO. We introduce a Hamiltonian of such 



t-V model as, 
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Here, Cj denotes the annihilation operator of fermions 
and rij(= c'-Cj) is its number operator. The index (ij) are 
the nearest neighbor(nn) pair sites. We deal with the tri- 
angular lattice with anisotropy in one of three directions. 
The geometry is reflected in t,t' and V, V as shown in 
FigEJ where the dashed values are the anisotropic ones. 

We first take the classical limit, t = t' = 0, for sim- 
plicity, which exactly corresponds to the Ising 5=1/2 
spin system. According to Wannier, the ground state 
of the regular triangle, V = V, is macroscopically de- 
generate and disordered, which are classified into two 
groups^. One group has the staggard alignment of 
charges (0101- • •) in one direction which stack to form 
regular as well as random stripes shown in Fig. ^ Here, 
and 1 denote the absence and the presence of charges, 
respectively. The other states are basically made up 
of three sublattices as shown in Fig. EJa). Two out of 
three sublattices(A and B) have the fixed population(l 
and 0) while each third sublattice site (C) can be arbi- 
trary. Thus the system is piled up by the hexagon unit 
in Fig. Efb), whose binding energy is always —3V. Be- 
sides the regular three sublattices, there are numbers of 
contingent freedom shown in Fig.|5fc). Since the energy 
of the three sublattice states including these contingents 
does not depend on the filling within 1/3 < p < 2/3, the 
system is gapless at p = 0.5[l^j. This is in sharp contrast 
to to the striped states with the excitation gap of 2V. 

When the anisotropy is introduced, the macroscopic 
degeneracy is lifted. At V > V , the vertical stripe 
type in Fig. ^a) becomes a unique ground state. As 
for V < V, the charges align in a staggered manner 
along the anisotropic direction, e.g. as in Fig. [Jc). Since 
the each chain can choose either (0101) or (1010)-type 
of alignment, there remains a semi-macroscopic 2-^ -fold 
degeneracy, and the disorder survives in the s-direction. 
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FIG. 1: Schematic description of the typical CO on a 
anisotropic triangular lattice. The staggard alignment of 
charges in one particular direction(CO chain) with either 
(0101) or (1010)-type of configuration is present. Represen- 
tative regular stacking of CO chains are (a) vertical and (b) 
diagonal striped states. One example of the random align- 
ment is given in (c). 
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FIG. 2: (a) Three sublattice structure; two sublattices( J 4 
and B) have the and 1 configurations, respectively, while 
the third one(C) could individually chose either or 1. (b) 
Hexagonal unit of the three sublattice structure. The center 
site (C) is free, surrounded by the staggard edges, (c) Repre- 
setative contingents. When the three nearest C-sites happend 
to have the same configuration, the center site{A or B in big 
circle) can switch its number without the energy loss. 



The classical binding energy, Ec, of the anisotropic tri- 
angle under the fixed value of V + V = Kum is given 
in Fig. EUa). The energy line follows Ec oc VN/2 and 
-VN/2 atV<V and V > V, respectively, and equals 
exactly with the three sublattice ones at V = V . 

In the highly frustrated disordered region an introduc- 
tion of small Uj is expected to lift the degeneracy. Let us 
start from the regular triangle at V = V ^> t = t'. The 
first term of eq.Q is considered as perturbations. In the 
striped states, the first order correction is absent and the 
energy changes to ~ (V — t 2 /V)N/2. Since this effect is 
small, the stripes are quite robust. 

In contrast, the three sublattice states are modified at 
the first order level. For demonstration, we first deal with 
the simplest p — 1/3 case. The charges align equally 
spaced like a Wigner crystal, e.g. in the A-sublattice. 
Next we add a single charge, and then there are 2N/3 
basis with the extra charge on one of the B or C-sites. 
The charges can hop by t to its neighboring B or C-sites 
as in Fig.^Ja) without the loss of binding energy. Resul- 



tantly, the charge can move around avoiding the A-site 
just like a pinball, so we call this "a pinball liquid" . The 
same situation holds at p = 2/3 where we regard the 
"pins" as a hole-rich site after the particle-hole transfor- 
mation. In the dilute limit, the extra charges thus form 
a honeycomb band. Such simple understanding is, how- 
ever, non-tribial when the charge density becomes close 
to half-filling, p ~ 0.5, where about one-half of the B- and 
C-sites are filled with charges. The "balls" can at least 
move from one site to another following the rule shown 
in Fig. 0{c) . Still, there are several issues to be exam- 
ined. First problem is that a significant number of states 
including the contingents in Fig. |3[c) will prevail, which 
might mix by t with regular three sublattice states and 
destroy the " pins" . Another thing is that the free pin- 
ball picture no longer holds since there are interactions 
between balls themselves. Also, the validity of the above 
discussion is not clear when the anisotropy is introduced. 

To clarify these points we proceed the numerical calcu- 
lations on eq.Q) at half-filling. We execute the exact di- 
agonalization at T=0 on A^=4x6=24 cluster. Then both 
the stripe type of two-fold states and the three sublattice 
states are compatible. Figure shows the structural fac- 
tor of the charge-charge correlation function at V=V' =10 
and t=t'=l under the periodic boundary condition. It is 
defined as, C k = ^ £ im ((n ; - l/2)(n m - l^))^ 1 ^. 
The three sublattice type of peak structures at k = 
(±7r, ±27r/3), (0, ±47r/3) are observed, which grow with 
increasing V. To see the details, we calculate the three- 
body correlation function, Pa(j) = {n\(l — ri2)nj), which 
reveals the population of the j-th site when the 1st 
and the 2-nd site is present and absent, respectively. 
The ones along the y-direction for several choices of 

V = V are shown in Fig. HJa). A clear three fold 
structure of A-B-C type with different density, namely 
Pz{A) +P 3 (C) ~ 2P 3 (B), is present, suggesting the ex- 
istence of particle- hole symmetry. Here, P3(A) denotes 
the correlation when j-th site belongs to the A-sublattice. 
Thus, although the amplitude of the "pins" are squeezed 
to ~ 0.8 due to the quantum effect, a characteristic 
"pinball" -type of three body correlation is present. It 
originates from the geometrical frustration and essen- 
tially differs from the two-body ones in the CO states. 
The correlation decays with distance so that we cannot 
determine whether it lasts in the bulk limit or not within 
our scheme. Nevertheless, the amplitude of "pins" are 
growing with V to form a firm structure at least over 
considerable distances. 

The pinballs sustain even when the anisotropy of V ^ 

V or t ^ t' is introduced. Actually, in the case of p = 1/3 
and 2/3 with a single "ball", the only modification is 
the binding energy per "ball" from the previous — 3V to 
— 2V — V. If only the first order correction of energy by 
t, Ei, is finite the three sublattice state should overcome 
the others as in Fig. E2a) and form a new phase at V ~ 
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FIG. 3: (a) Classical binding energy, Ec, at the fixed value 
of V + V' = Vsxaa as a function of V . The vertical and dis- 
ordered striped states are realized at V < V' and V > V' , 
respectively, while at V = V', the three sublattice states with 
Ec = (2V + V')N/6 becomes degenerate with these stripes. 
Energy correction, E\, up to first order perturbation of t low- 
ers the energy of the three sublattice state, (b) Ground state 
energy Eq of eq.JIJ at N=24 with V BU m = 20 and t = t' = 1. 
The diagram is classified into three regions; (I)vertical stripe, 
(Il)pinball liquid, and (Ill)disordered stripe, where Eq be- 
haves almost linear to respective Ec's. 



FIG. 4: (a) Representative basis of the three sublattice states 
near p = 1/3 that mixes in the first order of t at p = 1/3 
with one extra particle, (b) Schematic illustration of a pinball 
model, (c) Representative basis away from p — 1/3, that is 
allowed to mix by t. The "ball" can move from i-th to j-th 
site when the extra two neighboring sites of i-th site with big 
circles together have the same number of balls with those of 
the j-th site. 




V . Figure Gib) shows the ground state energy of eq.Q) 
together with Eq at fixed V + V = 20. The diagram 
is separated into three regions; the energy line behaves 
parallel to different classical states, (I) vertical stripe, (II) 
three sublattice, and (III) disordered stripe. The pinball 
liquid realized in (II) extends at V ^ V over a width of 
order t. Actually, the first order energy gain Ec — Eq 
at V — V is several times larger than the second order 
gain in the (I) and (III) regions. In this way, the pinball 
liquid is quite dominant at around V ~ V' . 

It is natural to expect that the gapless situation found 
in the classical three-sublattice states still holds when t 
is introduced. Since the direct gap estimation under the 
finite system size is difficult and unreliable, we instead fo- 
cus on the coherence length of each characteristic state. 
We twist the boundary condition as e 1 ^, where (f> = and 
7r correspond to the periodic and anti-periodic boundary 
conditions, respectively. The ground state energy, E(cj)), 
should have considerable (^-dependence only when the co- 
herence lasts longer than the system length. Figure Efa) 
shows the variation of E((f>)— E(0) as a function of 4> when 
twisting the boundary in the y-direction with V + V 
fixed. In the regions (I) and (III), E(<f) remains almost 
constant, reflecting the insulating character. In contrast, 
the ones in (II) show large <fi dependence, indicating that 
pinball wave function has a delocalized character at least 
over dozens of sites. Figure [7Jb) shows the amplitude of 
the energy variation, E(t:) — ^(O), when twisting <fi in 
the x and y-directions. One actually finds a significant 
difference that separates (II) from other two states. 

In summary, we proposed the existence of a pinball 
liquid state on a triangular lattice in the spinless fermion 



FIG. 5: (a) Charge structural factor, Cu, at 

V=V'=W,t=t'=l. (b) V(= V")-dependence of the peak at 
k= (0, ±4tt/3), (±7t, 2tt/3). 

systems with nn Coulomb interactions over a wide range 
of filling, l/3<p<2/3. It is an ordered liquid formed 
out of disorder by the lifting of macrscopic degeneracy 
and has a nobel three-body correlation. Thus, although 
the state is metallic it is driven by the local correlation 
of electrons at V ^S> t and is not related to the details 
of the Fermi surface, namely the geometrical structure of 
t. It is also insensitive to the anisotropy of V ^ V and 
survives quite away from the regular triangle. Further, it 
is expected to sustain when the degrees of freedom other 
than charges such as spins and orbitals are introduced. 

This nobel liquid is quite universal. For example in 
the hard core boson systems (XXZ spin systems), Py,{j) 
also shows a clear three-fold pinball character as given 
in Fig. Elb) . Therefore we conclude that it is not much 
influenced by the statistics and is also applicable to the 
large-JJ region of EHM at quarter-filling. 

Recently there had been other theoretical studies on 
the anisotropic triangular lattice that argues the exis- 
tence of the three- fold states E3 ■ Unfortunately, their 
treatment based on the mean-field wave functions rather 
take account of the instability of the Fermi liquid to- 
wards the two-body correlations, and could not capture 
the essential features of the pinball liquid, namely the 
frustration induced three-body correlations. The present 
theory gives the explicit and clear-cut answers to the no- 
bel nature of this system alluded to preliminarily in such 
previous studies |lflll8lll9| . 
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FIG. 6: Three-body charge correlation function, Ps(j), along 
the second chain(i = 7 — 12) at t = t' = 1 for several choices 
of V = V'. Here, i denotes the location of the site indicated 
in the inset. Those of the (a) the spinless fermion and (b) the 
hard core boson systems are shown. 
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FIG. 7: (a)Variation of the ground state energy, E{(f>) — 
E(0), as a function of <j> = — n, under the twisted boundary 
condition in the {/-direction for several choices of V with V + 
V' — 10 fixed, (b) The amplitude of variation E(-k)—E(0) as a 
function of V along V+V' = 10, when the boundary is twisted 
in the x- and y-directions. Regions (I)- (III) correspond to 
those of FigEH 



Related experimental studies are explored on the CO 
materials, 0-ET2X. Recently, a coexistence of the short 
range diffuse spot of two-fold and three-fold periodic- 
ity is observed in the metallic states of materials con- 
taining the Rb0| and Cs0l i° ns J which have V/V ~ 
0.85-0.9 20J . From our strong coupling viewpoint, both 
two-fold striped and pinball liquid states are compati- 
ble (though not coexisting) without having special as- 
sumptions such as long range interactions or the partic- 
ular shape of the Fermi surface. If the systems were in 
the critical region of these two competing characteris- 
tic states, some additional effects such as impurity, tem- 
perature(entropy effect), or electron- phonon interaction 
might induce a coexistence or domain structure that cope 
with these anomalous experimental findings |ll|. 

In conclusion, we found a new type of liquid state 
in the charge ordered system under geometrical frustra- 
tion, which is possibly relevant at temperature range of 
J < T -C V. In the strong coupling limit, the frustra- 
tion among electronic interactions, V, drives the system 
to disorder. The introduction of small but finite t(3> V) 



transforms the system to a strongly correlated quantum 
liquid with three sublattice correlation, A-B-C, whose 
electronic density is, + rig ~ 2n.o About half of the 
charges become a Wigner-crystal-solid while the rest re- 
mains a liquid. It is nobel in that an originally unique 
charge degrees of freedom spontaneously separates into 
the solid and liquid-like part and coexist in the same 
system. Such correlation induced order out of disorder 
distinctively differs from the usual charge ordering or 
charge-density-waves formed from the instability of the 
Fermi liquid. This state remains valid regardless of elec- 
tronic filling, the shape of the Fermi surface, anisotropy 
of V, and also the statistics of the particles, and endows 
a new detailed insight on the topic of frustration. 

We thank K. Kubo, S. Miyahara, I. Terasaki, R. Kondo 
and M. Watanabe for discussions. 
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